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The role of the auxihary scalar field (j> of Brans-Dicke theory played in baryon number asymmetry 
is discussed in this paper. We consider a derivative coupling of this gravitational scalar field to the 
baryon current b or the current of the baryon number minus lepton number b-l based on a 
series of works of R. Morganstern about the Brans-Dicke cosmology. We find that the spontaneous 
baryogenesis by this coupling is capable to yield a sufficient baryon asymmetry ns/s ~ 10^" for 
• the time of the grand unification is in a little advanced. In addition, Davoudiasl et al have recently 

' introduced a new type of interaction between the Ricci scalar R and the bayon current , d^RJ''' 

, and also proposed a mechanism for baryogenesis, the gravitational baryogenesis. However, the Ein- 

stein equation tell us that _R = in the radiation-dominated epoch of the standard FRW cosmology. 
In this paper we reconsider the feasibility of having gravitational baryongenesis with such a form of 
^ , interaction in radiation-filled Brans-Dicke cosmology. We will show that R does not vanish in this 

■ case and the required baryon number asymmetry can also be achieved. 
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I. INTRODUCTION 
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' The origin of the baryon number asymmetry remains a big puzzle in cosmology and particle physics. Conven- 

' tionally, it is argued that this asymmetry is generated from an initial baryon symmetry phase dynamically as long 
1 as the following conditions are satisfied (l)baryon number non-conserving interactions; (2) C and CP violations; 
(3)out of thermal equilibrium. When the CPT is violated dynamically, however the baryon number asymmetry can 
be generated in thermal equilibrium 0- In connecting to dark energy, a class of models of spontaneous leptogenesis 
O ' 0, 0, l^are investigated recently by introducing a interaction between the dynamical dark energy scalars and the 
, ordinary matter. Specifically, a derivative coupling of the quintessence scalar field Q to the baryon or lepton current 
is under considering: 



An* - d^Qjr (1) 



One silent feature of this scenario for baryogenesis is that the present accelerating expansion and the generation of 
' the matter and antimatter asymmetry of our universe is described in a unified way. 

Recently, a new mechanism called gravitational baryogenesis in thermal equilibrium has been proposed by 
Davoudiasl et al [^. They introduced explicitly an interaction between the Ricci scalar curvature with derivative 
and the baryon number current: 

^ = -^^^RJ'- (2) 



And the baryon number asymmetry is given by 

ub R 



s APT' 



(3) 



which shows that ns/s is determined by the value of R, however the Einstein equation, R = SttGT^^ = 8TrG{l — 3uj)p, 
tells us that R = in the radiation-dominated epoch of the standard Friedmann- Robertson- Walker (FRW) cosmology. 
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Davoudiasl et al in Ref.jH have considered three different possibilities of obtaining a non- vanishing R which include 
the effects of trace anomaly, reheating and introducing a non-thermal component with w > 1/3 dominant in the early 
universe. In the brane- world scenario Shiromizu and Koyama in Ref.Q provided another example for R example for 
R ^ 0. A generalized form of the derivative coupling of the Ricci scalar to the ordinary matter 

A„t ~ d^f{R)J^. (4) 

is also discussed in 3- Taking f{R) ^ InR, they have shown that df_if{R) ^ df^R/R does not vanish and the required 
baryon asymmetry can also be generated in the early universe. 

As it is well known, Brans-Dicke theory of gravity is one of the simplest extensions/modifications of Einstein's general 
relativitvlol ITol Hlf . Compared with general relativity, as well as the metric tensor of space-time which describes the 
geometry there is an auxiliary scalar field (f) which also describes the gravity. The role of the gravitational scalar 
in Brans-Dicke theory has been much exploited to achieve inflation, quintessence |1^ il^ JJ„ JJ,, 2&, 12/ 2^ 
and k-essence The main constraints on this gravitational scalar come from the solar system experiment |22j . In 
addition, the testing of Brans-Dicke theory using stellar distances, the CMB temperature and polarization anisotropy 
have also been discussed in [2^l23 |. In this regard, here we focus particularly on the scalar of the Brans-Dicke theory 
to model baryo(lepto)genesis behavior by combining both the spontaneous symmetry breaking and its gravitational 
origin in Brans-Dicke cosmology. 



II. SPONTANEOUS LEPTOGENESIS 



We consider a derivative coupling of the auxiliary gravitational scalar cf) to the ordinary matter. So we have now 
an effective Lagrangian 

r-eff = ^a,0J^ (5) 

where M is the cutoff scale which, for example, could be Planck mass Mpi or the scale of grand unification theory 
Mgut, and c is the coupling constant which characterizes the strength of gravitational scalar interacting with ordinary 
matter in the standard model of the electroweak theory. Taking = Jg, during the evolution of the spatial flat 
FRW universe, Ceff in the equation (5) gives rise to an effective chemical potential fib for baryons: 

c c • c • 

-j^d^(j)J^' ^ ^(l)nB^^(p{nb-ni), (6) 

= = -^fc. (7) 

In thermal equilibrium the baryon number asymmetry is given by (when T ^ rub) 



^b = ^[Y + 0{-)]^c-^, (8) 



where gb counts the internal degrees of freedom of the baryon. Using the familiar expression for entropy density 

« = '-^9*T\ (9) 



we arrive at the final expression for the baryon to entropy ratio: 

UB 15c gb4> 



s 47r2 g^APT' 

here 4> in Eq.(lO) can be obtained by solving the equation of motion of auxiliary scalar given below 



(10) 



(2w - 3)(^ -I- 3i7^) + = -^[UB + (0 + iH)nBl (11) 



where H is the Hubble constant and Tm is the energy-momentum scalar of the ordinary matter in the universe. 

In applying the Brans-Dicke theory to cosmology, we may write down the flat Robertson- Walker line element as 

ds^ = -dt^ + a'^{t)gijdx'dx\ (12) 
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Where i, j run from 1 to 3, a{t) is the scale of the non-compact 3-dimensional flat space. The action of Brans-Dicke 
theory reads 

+ J d^X^/^-hm, (13) 



here uj is the parameter of Brans-Dicke theory; Then the corresponding field equation is given by 

1 

2^ 



The field equation for (p reads 



'3 + 2uj 

Therefore, the fundamental equations of Brans-Dicke cosmology may be taken as 



^0;^'/';. + ^^5m-V,0V"0 (14) 



= -^T^- (15) 





St: p 




3 ^(j) ' 


d 


Stt (-3 


a 


3 (- 




Stt 






a 


-3 -1- 2uj 



2+^7; (17) 



-3 + 2tj)0 3 , 

(l-3u;)p. (18) 

Here w is the parameter of the state equation of the ordinary matter(?i; = p/ p). The exact solutions to Brans-Dicke 
cosmologies in flat Friedmann Universes had been drawn by R. Morganstern from the early 70's "2^, "2^, '23| • He found 
that as a — > the expansion parameter a has the same t^^^ behavior for all curvature k in radiation-filled Brans-Dicke 
cosmology 26] . Moreover, as t ^ it is found that the flat-space solutions display the same time dependence (and 
same dependence upon units) for a and (f> as do these same quantities in the radiation cases (fc = 0,±1) 27]. Therefore, 
we now turn to investigate the feasibility of having spontaneous baryogenesis in the early universe with Brans-Dicke 
gravity in the case of w = 1/3, so we have 

p oc a^^ and <p cx a^^. (19) 

In the research of Brans-Dicke cosmology, a power law form for both the scale-factor and scalar field has ever 
been applied to investigate the dynamics of a self-interacting Brans-Dicke field to account for the acceleration of the 
universe at late epochs t ^ to in |2^. More importantly, the exact general solutions to Brans-Dicke cosmology for 
p — p/3 found in |2y| has also demonstrated a limiting power-law behavior for early epochs t — > 0. For these reasons, 
we also consider a power-law form for the scale factor as an approximation to find its asymptotic behavior in the early 
universe, 

a = a^„(T^)". (20) 

Thus we have H = n/t, then Eq.(18) of w = 1/3 immediately yields the power law form for the scalar field 
(t> — 4>in{t/tin)^'^^'' and (j> = <pin(,t/tin)^~^" ' i_^3„ ■ Giving this result, now we can write down the parameter- 
ized Friedmann equation of (16) in a radiation-dominant epoch, 

^ STTjl -3n)p ^Un^,_,^ _ _ ^21) 

On the other hand, the radiation is in equilibrium so that p — ^g^T^ can be available in the final calculation. 
Therefore, when the interaction of the spontaneous baryogenesis is decoupled, a simple relation between the age and 
temperature of the Universe is obtained, 

t-'-'- = -. MljzM (i^„)-3«!l! (22) 

30„[n2 + n(l-3n)+c^(l-3n)2/6] ^ ^ 30^ 
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Substituting Eq.(22) into the power law form for the scalar field, we obtain 



(t'i7iit/tin) 



-3n 



''^^'-'"^ 49.T\ (23) 



3[n2 + n(l - 3n) +w(l - 3n)2/6] 30' 



As far as the radiation-filled Brans-Dicke cosmology is concerned, R. Morganstern has found the limiting behavior 
of the exact solution for the expansion parameter is oft^^^. Considered spontaneous baryogenesis takes place at the 
early times of the universe, when Eq.(21) possesses an asymptotic solution as n — > 1/3, so it requires that 

-uj{l~3nf/6^n'^+n{l-3n). (24) 

in the limiting behavior. The equal sign of above equation would come into existence only when n = 1/3 in our model. 
However, to say at least, Eq.(24) is still a fine approximation when we consider the exact solution to radiation-filled 
Brans-Dicke cosmology which possesses an asymptotic behavior of n ^ 1/3 at the early epoches. Therefore, the ratio 
of baryon number to entropy in this case is given by 

ns ^ 15c SiT^gi, (1 - 3n) 

s ~ 47r2 ' 90Af2 ' n2 + n(l - 3n) + w(l - 3n)2/6 ' ' ^ ' 

where t and T refer to the age and temperature of the Universe, respectively. 

In the scenario of spontaneous baryogenesis, the baryon asymmetry is generated in thermal equilibrium. This 
requires that baryon number violating interactions occur rapidly (F^ > H). However, if the i?- violating interactions 
keep in equilibrium until 0-^-0, the final baryon asymmetry will be zero. Denoting the epoch when the i?-violating 
interactions freeze out by (corresponding to to), i.e., Tsitij) = H{To), the final baryon number asymmetry is 
obtained 

^ _ -, (1 - 3^^) tpTp^ 



-n(l-3n)-f w(l-3n)2/6 ' NP 



In the above numerical calculations, we have used gi, ^ C'(l)- Furthermore, we may a priori assume that 
ni+n(i-^ri)+^{i-3n)'i /6 ~ C(10~^) for Considering \uj\ > 500. However, as we shall see in the last of this section, 
it will not change our result much if we assume this expression takes the value in the order ~ 10^"^. After the in- 
teraction of the spontaneous baryogenesis decouples the universe will continue to evolve into a radiation-dominated 
epoch until the primordial nucleosynthesis takes place. But we must note that thereafter the evolutional behavior of 
the scale factor will deviate from the n ^ 1/3 law along with the increasing of the time t. Therefore, we must be 
interested in the epoches of the back direction from this decoupling time. As a very good approximation, we can recall 
the relation of equilibrium thermodynamics that ^ = (^)4/4n _ ^p^-^-j:^ _ (£^£i/|0^)-^ obtain the evolution 
of the temperature. To make our discussion to be as independent as possible of the cosmological evolution models, 
here we choose the values of Tqut at the grand unification phaser transitions epoch to be the input parameters in 
our model. Firstly, we choose Tqut ^ lO^^Gew considered from particle physics. Then the induced baryon number 
asymmetry in radiation-filled Brans-Dicke cosmology is 



— 3- 



S A/2 

142^^,3 ^GUT 



*l/3 



.l..l„.W.^, (27) 

here we have taken the cut-off factor M ~ Mp/. Then yielding a sufficient ^ ~ 10"""^° would require that tcc/T ^ 
10~"'^^Gew~^ ~ 10~^^sec. This value is not very far away from the Planck time (t^; ^ 10~''^sec) but still can be 
achieved, just need the grand unification phase transition to take place in a little advanced comparing with that of 
Einstein's cosmology. If we consider that the index of the power law n will show a preference of a little departure from 
1 /3 along with the decrease of the temperature until Tp , then the requirement on the time when the grand unification 
becomes to take place would be further released to be > 10~'^*sec. Along this way, the spontaneous baryogenesis 
using the Brans-Dicke scalar is also possible to yield sufficient ub/ s at a radiation-dominated epoch. 

On the other hand, if the B-violating interactions conserve B — L, the asymmetry generated will be erased by the 
eletroweak Sphaleron [2^. Hence, now we turn to leptogenesis [s^,!^- We take in Eq.(5) to be J^^b-l- Doing 
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the calculations with the same procedure as above for = J'^s we have the final asymmetry of the baryon number 
minus lepton number: 

ns-L _ (1 - 3») tpTp^ 

s ~ ' n^ + n{l-3n)+uj{l-3ny/6' A'P ' ^ ' 

The asymmetry ub-l in (28) will be converted to baryon number asymmetry when electroweak Sphaleron B + L 
interaction is in thermal equilibrium which happens for temperature in the range of lO^GeV ~ lO^^GeV . Tp in (28) 
is the temperature below which the B — L interactions freeze out. 

In the Standard Model of the electroweak theory, B — L symmetry is exactly conserved, however many models 
beyond the standard model, such as Left-Right symmetric model predict the violation of the B — L symmetry. In this 
paper we take an effective Lagrangian approach and parameterize the B — L violation by higher dimensional operators. 
There are many operators which violate B ^ L symmetry, however at dimension 5 there is only one operator, 

Cc = jIlIiXX + H.c, (29) 

where / is a scale of new physics beyond the Standard Model which generates the B — L violations, II and x ^I'e the 
left-handed lepton and Higgs doublets respectively, when the Higgs field gets a vacuum expectation value < x >^ v, 

2 

the left-handed neutrino receives a majorana mass ~ 

In the early universe the lepton number violating rate induced by the interaction in (29) is jH3 



■I! 



Tl - 0.04—. (30) 



Since is proportional to T^, for a given /, namely the neutrino mass, B — L violation will be more efficient at 
high temperature than at low temperature. Requiring this rate be larger than the universe expansion rate ^ n/t ^ 

z[n-+n{i-ln]+lt-^n)-/6\ ' i^9*T^ ^^^^ temperature To, we have 

/ - i — tDT^D)^\n-.l/3 - {O.UtDT^D)^. (31) 

n 

Therefore, we obtain a To-dependent lower limit on the neutrino mass: 

6 X W^Gev^ 



- lev, (32) 

where we have using the energy scale of electroweak SSB phase transition of about 2A6Gev. Taking three neutrino 
masses to be approximately degenerated, i.e., mi ~ m2 ~ ma ~ rh and defining ^ = 3rh, one can see that 
for toT^D ~ 3 X lO'^^Gev^ GeV, three neutrinos are expected to have masses m around 0{leV). The current 
cosmological limit comes from WMAP [s^ and SDSS The analysis of Ref.fssj gives E < 0.69eV. The analysis 
from SDSS shows, however that S < 1.7eV [33|. These limits on the neutrino masses requires toT^D be larger than 
5.7 X IQ^^Gev^ and 9.2 x lO^^Gew^. The almost degenerate neutrino masses required by the leptogenesis of this model 
will induce a rate of the neutrinoless double beta decays accessible for the experimental sensitivity in the near future 
135]. Interestingly, a recent study ,36] on the cosmological data showed a preference for neutrinos with degenerate 
masses in this range. Therefore, we may recall the final expression (26) for spontaneous baryogenesis by using the 
Brans-Dicke scalar. If the cut-off factor M takes the value from the grand unification scale 10^^ GeV to the Planck scale 
10^^ GeV, then yielding a sufficient baryon asymmetry may requires that ^ x 10~^ < ^'^ j^n(i~l^+u{i-3n)2 /e, ^ '^l'^- 
It is quite easy to be achieved for a large number of Brans-Dicke parameter uj. 

The experimental CPT test with a spin-polarized torsion pendulum |37j | puts strong limits on the axial vector 
background h^^ defined by £ = b^e'^'^^f^e 

\b\ < IQ-^^GeV. (33) 

For the time component bo, the bound is relaxed to be at the level of 10~^^ GeV js^. In our model, assuming the 
auxiliary scalar couples to the electron axial current the same as Eq.(5), we can estimate the CPT- violation effect on 
the laboratory experiments. As we know, the inverse of the Newtonian gravitational constant G takes the average 



6 



value of the scalar 0, in the spirit of Brans-Dicke theory. The variation of G in the present time is bounded as 
g < 1.2 X IQ-'^^Gev by 113, 113. Hence the induced CPT- violating bo is 

< -1.2 X lO^^^Gey, (34) 
which is much below the current experimental limits. 



III. GRAVITATIONAL LEPTOGENESIS 

In this section we will recall the gravitational baryogenesis 6] in thermal equilibrium, which may fails to work in 
the framework of Einstein's cosmology just has been introduced in the beginning of this paper. So we investigate it 
again in Brans-Dicke cosmology. The effective Lagrangian of a CP- violating interaction between the derivative of the 
Ricci scalar curvature R and the baryon number(B) current J'^ still have the form of Eq.(2). Taking = J'^, during 
the evolution of the spatial flat FRW universe, Ceff in the equation (2) gives rise to an effective chemical potential 
fib for baryons: 

-^^^RJ' ^ -M^^^^ - -M^^i^" ~ ^i)^ (35) 

= c-^ = -fxi. (36) 
We could directly give out the final expression for the baryon to entropy ratio in thermal equilibrium: 

ub 15c gbR 



R in Eq. (37) can be obtained by solving the equation of motion of the Ricci scalar given below 



(37) 



^ - - + 2c.(^ - ^) - 3(- - ^) - 9Hi- - ^) - 9H- (38) 

= 6^ + 6^^-12^. (39) 

Obviously, there are two methods which may be used to find the behavior of R of the asymptotic solution {n — s- 1/3). 
The first simply involves substituting the power law approximation of (20) in the function of the curvature scalar (39) 
of the Robertson- Walker metric, 

R = I2n[{n - if - n^] ■ t-^ . (40) 

The second method can be used even when the explicit solutions are not known, and moreover it provides a convenient 
check on the asymptotic behavior of the exact solution (namely the assumed power law form (20) when n — s- 1/3). 
By substituting the power law of the scale factor in Eq. (38) (which is resulted directly from the field equation in 
Brans-Dicke cosmology (14)), we again have 

R = -2(1 - 3nfLj ■ r^. (41) 

On the other hand, here we should recall the relation (24) as a fine approximation at the earliest times. After that, 
we can compute the limiting behavior of R using both methods, and of course wc find that the results agree. Taking 
n — > 1/3, we immediately have R = 12f^^/9. Thus we obtain the ratio of baryon number to entropy, which is given 
by 

ub 15c 12(7b 



s 47r2 %g^NPT 



■ (42) 



where t and T also refer to the age and temperature of the Universe. A rapidly violating of baryon number is still 
required in this case. In addition, if the _B-violating interactions keep in equilibrium until i? — > 0, the final baryon 
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asymmetry will be zero. Denoting the epoch when the _B- violating interactions freeze out by T^i (corresponding to 
tu), i.e., Tsito) = H{T£)), the final baryon number asymmetry in this case may be given by 

UB 15c I2gb ^ _3 
— ~ ■ ■ to 



s 47r2 ' 9g^APTD 



In the numerical calculations above, we have also used gb ^ and g^ ~ 0(100). Recalling the equilibrium 

thermodynamics tu ^ tGUT{'^-^^)~ at a radiation-dominated epoch w — 1/3, then the final result will be given 
under the consideration of the limiting behavior of the exact solution in , 

ns c ^^_2 -3 {W^Gev)-^TD~ 



X 10 -t GUT l^rPFZ l"^i/3 



-2 ^-3 

2 ' ""^ M'^Td 

C i^-137^„„,-9 .-3 . 



^ X 10- Gev~^ -t-^GUT-^. (44) 

Similarly, if the B-violating interactions conserve B ~ L, the asymmetry generated by gravitational CP-violation will 
also be erased by the eletroweak Sphaleron 22] . So we again have to turn to leptogenesis js^, |^ ■ We take J'' in 
Eq.(35) to be J^b-l- Doing the calculations with the same procedure as above for J^^ = J^b we have the final 
asymmetry of the baryon number minus lepton number in a similar way to spontaneous leptogensis 

^ ^ X lO-'^'Gev-' ■ r'avT ■ (45) 

Taking c ~ 0(1), ub-l/s ^ 10"^° and M ~ lO^^GeF, for the desired result is obtained when To < lO^^Gev, then 
the time tcuT when the grand unification becomes to take place may be required to be 

tauT < IQ-^^sec. (46) 

That is to say, the gravitational baryogenesis in Brans-Dicke cosmology also requires the grand unification phase 
transition take place in a little advanced comparing with that of Einstein's cosmology. 

In the same way, we investigate the CPT-violation effect on the laboratory experiments for the coupling between 
the curvature scalar and the electron axial current. Firstly, viewed from Taylor expansion, when a variable quantity 
evolves to a great number, the derivative of a function of this quantity would tend to be more and more smooth 
and flatted along with the increasing of the rank of the derivative. Thus we have <C <C at the present 
time(here the infinitesimal unit At ~ Isec is chosen to obtain a large number for the age of the present universe 
relative to considered time interval). Secondly, the Bubble's parameter at the present time has been given as Hq ~ 
2 X IQ-'^^Gev. Therefore the induced CPT-violating bo in this case may be estimated at 



c • c , aa a , a 

t^to 



bo - —^Ro-jj^[-6^ + 6-~l2H' + l2--H]\t^ 



< -^\2A X lO^^^^Gei;^ • sec^^ + 1.2 x IQ-^^Gev ■ sec^^ - 9.6 x IQ-^^^Gev^] 

< 4.8 X 10-i22(^gy^ (47) 
which also is much below the current experimental limits. 



IV. CONCLUSION 



In summary, we propose in this paper a scenario of spontaneous baryogenesis in Brans-Dicke cosmology by in- 
troducing a derivative coupling of the auxiliary scalar cf) to the ordinary matter. Our model is capable to explain 
the baryon number asymmetry ub/s ^ 10^^" without conflicting with experimental tests on CPT. As a complemen- 
tarity, we also investigate the scenario of gravitational baryogenesis in radiation-flUed Brans-Dicke cosmology from 
a derivative coupling of the Ricci scalar to the ordinary matter introduced by Davoudiasl et al in this case the 
current baryon asymmetry can also be achieved at early epoches just for a non-vanishing R, in different to the model 
in Einstein's cosmology. As a result, we should also keep in mind, to obtain a sufficient baryon asymmetry in both 
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these two models in Brans-Dicke cosmology, the time when the grand unification becomes to take place would be 
required to be in a little advanced comparing with that of Einstein's cosmology. 

After this work was firstly submitted a related independent workji^l appeared in the archives, which has similar 
motivations^. 
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